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Abstract 

A fundamental problem in the theory of PT-invariant quantum systems is to determine 
whether a given system "respects" this symmetry or not. If not, the system usually develops 
non-real eigenvalues. It is shown in this contribution how to algorithmically detect the 
existence of complex eigenvalues for a given PT-symmetric matrix. The procedure uses 
classical results from stability theory which qualitatively locate the zeros of real polynomials 
in the complex plane. The interest and value of the present approach lies in the fact that it 
avoids diagonalization of the Hamiltonian at hand. 

1 Motivation 

When dealing with a non-hermitean operator such as 

H=p^ + ix^, (1) 
one needs to address two questions which do not arise for a Hermitean operator: 
Ql: Is the operator H is diagonalizable? 
Q2: Does the operator H have real eigenvalues only? 

For a randomly picked non-hermitean operator the answers to both questions are unlikely to be 
positive: it will neither have a complete set of eigenfunctions nor a real spectrum. However, 
operators with PT-symmetry pQ, 

[H,PT] = 0, (2) 

invariant under simultaneous application of parity P and time-reversal T, behave somewhat 'bet- 
ter.' PT-invariant operators tend to be diagonalizable but for the rare occurrence of exceptional 
points, and each of their eigenvalues must be either real or have a complex conjugate counter- 
part. Positive answers to Ql and Q2 are necessary in order to attempt a consistent quantum 
mechanical interpretation of the operator H since it can be similar to a hermitean operator only 
then [g. 

To answer these question for a given PT-invariant Hamiltonian is by no means straightforward. 
It is known, for example, that the operator H in Eq. JQ) does have only real eigenvalues W while 
the (likely) completeness of its eigenfunctions has, apparently, not yet been established rigorously. 
Perturbative results allow one to confirm that the spectrum of an initially hermitean operator such 
as the Hamiltonian of a particle in an oscillator-type potential remains real if a sufficiently weak 
PT-symmetric term is added 0]. As long as no degeneracies develop, this approach also makes 
plausible the existence of a complete set of eigenfunctions; they are, however, not necessarily 
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pairwise orthogonal with respect to the standard scalar product in Hilbert space. Technically, 
the difficulties are due to the fact that the cubic term in Eq. is unbounded on the real 
axis, and the unperturbed operator does not provide a bound for it. When restricted to a finite 
interval, a perturbation such as igx^,g £ R, is bounded, and one can reach general conclusions 
when perturbing the hermitean boundary value problem with a non-hermitean PT-symmetric 
term. Upon treating the PT-symmetrically perturbed square-well potential [S] in a Krein space 
setting, one can show |^ that its eigenvalues remain real if the perturbation does not move the 
(non-degenerate) real eigenvalues far enough along the real axis to create a degeneracy, which is 
necessary for complex eigenvalues to emerge. A similar result also follows by a non-perturbative 
approach when a "slightly" non-selfadjoint term is added to a self-adjoint operator as is described 
in 13. 

More is known for PT-symmetric systems with a /inite-dimensional state space which are 
described by complex symmetric matrices M. Let us consider an example which exhibits the 
essential features: the discretized PT-symmetric square well .8^. This model, sketched in Fig. ^ 
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Figure 1: Discretized PT-symmetric well: the wave function takes non-zero values at three points 
a; = 0, ±L only (cf. text) 



is obtained upon discretizing the configuration space of a particle moving freely between walls 
at a; = ±2L, subjected to a piecewise constant potential ±iZ, Z G R. Defining a wave function 
which takes values at the points a; = 0, ±L, ±2L, and satisfies "hard" boundary conditions at 
X = ±2L, an effective Hamiltonian is obtained. 



H = 1 1 \ , £, = 2mL^Z/h'^. (3) 




This matrix is invariant under the action of parity P, a matrix with ones along the minor diag- 
onal and zeros elsewhere, followed by complex conjugation, overall equivalent to Eq. (jJJ. The 
eigenvalues of H are given by the roots of its characteristic polynomial, 

PHi\) = \i\'~i2~e)), (4) 



reading explicitly. 



The possibility to analytically determine the eigenvalues of H provides immediate and exhaustive 
answers to both Ql and Q2, summarized briefly now. The zero eigenvalue (with its associated 
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eigenstate) persists for all values of ^, while the remaining two change their character with varying 
strength of the parameter Z. Three regions can be identified (cf. Part (b) of Fig. Q): depending 
on the magnitude of ^, there is either a pair of complex-conjugate or a pair of real eigenvalues. 
However, the matrix H is not diagonalizable for ^ — ±^/2: only a single eigenvector is associated 
with E± = 0, while the algebraic multiplicity of this eigenvalue is two 

For matrices M of larger dimensions no analytic expressions for the eigenvalues exist. To 
overcome this shortcoming, an algorithm has been proposed which is capable to detect whether a 
PT-invariant matrix is diagonalizable or not . The relevant information is coded in the minimal 
polynomial of the matrix which one can construct without knowing the eigenvalues of M, just as 
its characteristic polynomial. This approach answers Ql systematically, circumventing the need 
to numerically calculate the eigenvalues of M. This is important from a conceptual point of view. 

In the present contribution, a second, independent algorithm will be presented which answers 
Q2 for any PT-symmetric matrix. Both the number of its real eigenvalues and the number of 
pairs of complex eigenvalues are obtained by manipulating the coefficients of the characteristic 
polynomial of M. This information will be called the qualitative spectrum of M. 

The interest of the method proposed is due to the fact that it is possible to extract nothing 
but the desired information about the eigenvalues, namely their location relative to the real axis 
in the complex plane. Problems of this type arise when the stability of dynamical systems is 
addressed where it is crucial to determine whether the eigenvalues of a given matrix have negative 
real parts. 

In Section 13 the notion of inertia is introduced for Hermitean matrices, followed by Jacobi's 
criterion of stability for such matrices. Then, the breaking (or not) of PT-symmetry is described 
in terms of a modified inertia. Next, a theorem by Jacobi and Borhard is presented which locates 
the zeros of real polynomials in the complex plane. Section |3| combines all this to formulate 
an algorithm which, given a PT-invariant (or quasi-Hermitean) matrix, outputs the number of 
its real and of its complex eigenvalues. Finally, the algorithm is illustrated by applying it to 
the discretized PT-symmetric square-well potential introduced above, outputting correctly its 
qualitative spectrum. 



2 Stability and inertia of matrices 

Consider a dynamical system which is described exactly or, after some approximation, by the 
equation 

§ = M-x. (6, 

where M is a fixed hermitean (or real symmetric) matrix of dimensions (N x A'^) , and the vector 
x(t) gives the state of the system at time t. In many applications, one needs to know whether 
the solutions of Eq. ^ are stable: this is the case if all eigenvalues M„ of M have negative real 
parts, 

9leM„<0. (7) 

Indeed, no solution of Eq. © will grow without bounds if Q holds, making it possible to 
qualitatively predict the system's long-term behaviour. Let us characterize a matrix M by a 
triple of non-negative integers, its inertia 11 Qj with respect to the imaginary axis. 

In M = {v, S, n} , (8) 

where v and tt are the number of its eigenvalues with negative and positive real parts, respectively, 
while 5 counts the eigenvalues on the imaginary axis (cf. ^ for an illustration). A stable matrix 
M has an inertia of the form 

In M = {iV, 0, 0} , (9) 
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while a matrix is called marginally stable if none of its eigenvalues have a negative real part, 
allowing for the presence of purely imaginary eigenvalues, 

InM = {iV-m,m,0}, < m < iV . (10) 

Whenever tt > 0, the matrix M is called unstable since there is at least one solution of ^ which 
will grow without bound. 

2.1 Inertia of Hermitean matrices: Jacobi's method 

Jacobi devised an ingenious method PTT) to determine the inertia of a given (non-singular) her- 
mitean matrix L of size (TV x N). First, calculate the determinants d„ of its TV leading principal 
submatrices Li, L2, . . . , L^r = L, 

d„=detL„, n=l,...,iV, (11) 

all of which must be different from zero; second, write down a "+" followed by the sequence of 
signs cr„ of the N determinants d„. 



+ , fJi, (72, . . . , (Tat , cr„ = — ^=±1. (12) 

These (A'^+l) signs encode the inertia of the matrix L: the number of sign changes in this sequence 
equals the number tt of eigenvalues with positive real part, while the number of constancies in 
signs equals the number ly of its negative eigenvalues: 



# of constancies in (|12|l = tt 

# of alterations in (|12|l = i/ 



In L = (z^, 0, tt) . (13) 



The matrix L cannot have a zero eigenvalue, that is, 6=0, since all leading subdeterniinants 
including djv have been assumed to be nonzero. 

The following paragraph will show that it is possible to detect the location of the eigenvalues 
of a PT-symmetric (hence non-hermitean) matrix relative to the real axis by similar methods. 

2.2 Stability and inertia of PT-invariant matrices 

A non-hermitean matrix H with PT-symmctry satisfies ^ which implies that its characteristic 
polynomial 

N 

Ph(A) = ^/i„A" (14) 

n=0 

has real coefficients /i„ only, 

P^,(A)=ph(A*). (15) 

As a consequence, the zeros of this polynomial are either real or they come in complex-conjugate 
pairs. To distinguish between broken and unbroken PT-symmetry, it is useful to introduce the 
inertia of a matrix H with respect to the real axis, 

InsR H = {i/sR, (5sft, TTsff} , (16) 

where the triple {i'sr, <5sr, ttsr} of integers denotes the number of eigenvalues of H with negative, 
vanishing, and positive imaginary part (cf. Fig. lO)- The inertia of a matrix with real eigenvalues 
only, corresponding to unbroken PT-symmetry, reads 

InsiH = {0,7V,0}, (17) 
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Figure 2: Inertia of a (17 x 17) matrix with broken PT-symmetry; its imaginary and real inertia 
are given by InM — {5, 2, 10} and Inj} M = {6, 5, 6}, respectively. 



while broken PT-symmetry is signaled by an inertia of the form 

InjfH = {TO,7V-2m,TO}, m>0, (18) 

corresponding to m pairs of complex eigenvalues and {N — 2m) real ones. Let us know turn to 
the question how to determine the real inertia of a matrix with PT-symmetry. 



2.3 Zeros of real polynomials 



Given a real polynomial p{X) of degree N, Borhard ^2] and Jacobi JHI propose to proceed as 
follows to obtain the number of its real zeros. To begin, one determines the first (2iV — 2) Newton 
sums associated with the polynomial p{X) defined by 



so = iV, 



sn ^ X^l; + . . . + X 



N J 



n= 1,2,, 



, 2iV - 2 . 



(19) 



This is possible without knowing the zeros Ai, . . . , Aat, since one can JHl either define the numbers 
Sn recursively in terms of the coefficients hn of the polynomial or generate them by means of the 
identity 

dpiX) 



dX 



(soA ^ + siA ^ 



■MA). 



(20) 



Once the Newton sums have been calculated, one introduces the real symmetric (and Hermitean) 
matrix 

/ So Si S2 ■ ■ ■ Spf-l \ 

51 S2 ■■■ Sat 

52 Sn+1 . (21) 



\ sn-1 Sn 



S2JV-2 / 
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which is of Hankel type. One can thus apply the method presented in Section |2. II to determine 
its imaginary inertia. This is useful since Borhard ^2 and Jacobi 13 have shown^ that the 
inertia of 9)p, in fact, encodes the structure of the zeros of the polynomial p{X): 

,,,, { -K — V different real zeros, 
in io, = 0, ^ p(A) has | ^ ^.^^^^^^ ^^.^^ compiex-conjugate zeros. ^^^^ 

Let us imagine that the real polynomial p(A) is the characteristic polynomial ph{X) associated 
with a PT- invariant matrix H. Then, the result H22|l says that H has v pairs of different complex 
eigenvalues and (tt — v) different real eigenvalues if the Hankel matrix S)h associated with ph{\) 
has i^(7r) eigenvalues with negative (positive) real part. Expressed in terms of inertias, this result 
reads 

In f)H — {v,Q,it} iTi'g^U — — v,i/} . (23) 

The next Section will collect the results obtained so far and present them as an algorithm to 
determine the number of complex pairs and real eigenvalues of a PT-invariant matrix. 



3 Algorithm detecting complex eigenvalues 

Given a matrix H of dimensions {N x N) which is invariant under the combined action of parity 
P and time reversal T, Eq. (jJJ, here is an algorithm qualitatively determines its qualitative 
spectrum: 

1. Calculate the characteristic polynomial ph{\) of the matrix H; 

2. Determine the first {2N ~ 2) Newton sums s„ associated with the polynomial ph(A); 

3. Write down the Hankel matrix f)H l|21|l . defined in terms of the sums s„; 

4. Obtain the number of constancies tt and alterations v in the sequence of signs ()12|l giving 
the inertia of S^h as Ini^H = {i^, 0,7r}; 

5. Then, the inertia of the PT-invariant matrix H follows from the inertia Ini^n using H23|l 
with N = + 

Inji Y\^{v,N - 2v, v} . (24) 

Consequently, PT-symmetry is broken if > 0, and H will have v pairs of complex conjugate 
eigenvalues while the remaining [N — 2^) ones are real. Thu, the main result of this paper has 
been established. 



3.1 Example: The discretized PT-symmetric square well 

Let us work through the algorithm to detect the qualitative spectrum of the PT-symmetric 
discretized square-well potential described by the matrix H in Q — this time without solving for 
its eigenvalues. The derivative of its characteristic polynomial (0J reads 

%i^=3A2-(2-a- (25) 

Compare the expansion 



^H^-^) = 3A-1 + 2(2 - ^2)A-3 + 2(2 - ^'^f\-'^ + ©(A^^) (26) 



Ph(A) 



^The content of Refs. il21 and 1131 is described in IIUI . 
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Figure 3: Comparison of (a) the qualitative spectrum H31|) obtained algorithmically with (b) the 
exact eigenvalues of the Hamiltonian ^ describing the discretized PT-symmetric well 



with H20|l . to read off the first five Newton sums. The Hankel matrix associated with H is given 

by 

/ 3/2 {2-e) \ 

= 2 {2~e) , (27) 

\{2-e) {2-er J 

and its leading principal minors have determinants 

di=3, d2 = 6(2-e'), da = 20(2-^2)3^ (28) 

Depending on the value of the parameter ^, two different sequences of signs arise: for < 2, 
one has all dn positive, resulting in three constancies and no alteration: 

+ + ++ Ini^H = {0,0,3}, (29) 

while (^2 and turn negative for > 2, implying that 

+ + ^ Ini^H ={1,0,2}. (30) 

Using the relation (|24|l . the inertia of H with respect to the real axis is finally given by 

r {0,3,0} if iei<V2, 

'""^^i {1,1,1} if iei>V2. ^''^ 

Thus, the spectrum of H is real for < 2, while a pair of complex eigenvalues exists whenever 
> 0. This agrees with the exact result Q as depicted in Fig. 13 For ^ — ±2, the method 
cannot be applied since the matrix i3H in (|27|l develops leading principal minors with vanishing 
determinant. This is consistent with the fact that for these values of ^ the properties of the 
matrix H undergo qualitative changes such as the 'disappearance' of an eigenstate. However, this 
does not put the current approach in jeopardy since these exceptional points can be identified 
beforehand by running the algorithm presented in (^l which checks whether a given PT-invariant 
matrix is diagonalizable. In the present example, the points ^ = ±2 would be flagged as shown 
explicitly in 

In fact, modifications of the current approach have been developed by Gundelfinger and 
Frobenius (cf. ^U]) which are able to cope with the presence of at most three consecutive 
vanishing determinants d„. The general relation between vanishing principal sub-determinants 
dn of the Hankel matrix i^n and the zeros of the polynomial ph(A) is not obvious. In view of the 
example discussed above, it seems reasonable to conjecture that there is a close link between the 
non-diagonalizability of the matrix H and the existence of vanishing leading submatrices dn of 
i3H. 
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4 Discussion and Outlook 

An algorithm has been presented which is capable to determine whether the eigenvalues of a PT- 
invariant matrix H (or possibly a family of such matrices depending smoothly on parameters) are 
complex or not. It complements an earlier algorithm which detects whether a PT- invariant 
matrix does have a complete set of eigenstates. Thus, the fundamental questions Ql and Q2 
about PT-invariant systems spelled out in the introduction can be answered in a systematic way 
if the system is described by a matrix of finite but arbitrarily large dimension. 

Although desirable, it is not obvious how to generalize the algorithm presented here to op- 
erators such as H = + ix^ , acting in an infinite- dimensional space. This observation also 
applies to the algorithm for diagonalizability. From a numerical point of view, the algorithm is 
not particularly efhcient since a total of 2N determinants of order up to need to be calculated. 
However, the method proposed here is exact, contrary to any numerical implementation which 
would directly calculate the (approximate) eigenvalues of the matrix H. More efficient algorithms 
to determine the spectrum of a PT-invariant matrix are likely to exist — Sturmian sequences based 
on the Euclidean algorithm for polynomials |14j being the most promising candidates. 
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